Quantum State Tomography with Joint SIC POMs and Product SIC POMs 
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We introduce random matrix theory to study the tomographic efficiency of a wide class of mea- 
surements constructed out of weighted 2-designs, including symmetric informationally complete 
(SIC) probability operator measurements (POMs). In particular, we derive analytic formulae for 
the mean Hilbert- Schmidt distance and the mean trace distance between the estimator and the true 
state, which clearly show the difference between the scaling behaviors of the two error measures with 
the dimension of the Hilbert space. We then prove that the product SIC POMs — the multipartite 
analogue of the SIC POMs — are optimal among all product measurements in the same sense as the 
SIC POMs are optimal among all joint measurements. We further show that, for bipartite systems, 
there is only a marginal efficiency advantage of the joint SIC POMs over the product SIC POMs. In 
marked contrast, for multipartite systems, the efficiency advantage of the joint SIC POMs increases 
exponentially with the number of parties. 

PACS numbers: 03.67.-a, 03.65.Wj 



I. INTRODUCTION 

Quantum state tomography is a procedure for infer- 
ring the state of a quantum system from generalized 
measurements. It is a primitive of quantum computa- 
tion, quantum communication, and quantum cryptogra- 
phy, because all these tasks rely heavily on our ability 
to determine the state of a quantum system at various 
stages. One of the main challenges in quantum state 
tomography is to reconstruct a generic unknown quan- 
tum state as efficiently as possible and to determine the 
resources necessary to achieve a given accuracy, which 
can be quantified by various figures of merit, such as the 
trace distance, the Hilbert-Schmidt (HS) distance, or the 
fidelity [H,l. 

A generalized measurement in quantum mechanics is 
known as a probability operator measurement (POM). 
A measurement is informationally complete (IC) if any 
state is determined completely by the measurement 
statistics [3|-|5j]. In a d-dimensional Hilbert space, an IC 
measurement consists of at least d 2 outcomes, whereas 
a minimal IC measurement consists of no more than d 2 
outcomes. A particularly appealing choice of IC measure- 
ments are those constructed out of weighted 2- designs, 
called tight IC measurements according to Scott pi, 0] ■ 
Under linear quantum state tomography, they not only 
feature a simple state reconstruction formula but also 
minimize the mean squared error (MSE) — the mean 
square HS distance between the estimator and the true 
state @ . The construction of tight IC measurements has 
also been discussed in detail in Ref. Q. 

A symmetric informationally complete (SIC) POM 0- 
ITllJ is a very special tight IC measurement composed of 
d 2 subnormalized projectors onto pure states with equal 
pairwise fidelity of l/(d+ 1). They may be considered as 
fiducial measurements for state tomography for reasons 
of their high symmetry and high tomographic efficiency 



U2j . It is widely believed that SIC POMs exist in 
any Hilbert space of finite dimension since Zauner posed 
the conjecture although no rigorous proof is known. 
Analytical solutions of SIC POMs are known for d < 16 
and d = 19, 24, 28, 35, 48, see Refs. (MIDI and the ref- 
erences therein; numerical solutions with high precision 
have also been found up to d = 67 0, [111, [lj] ■ In addi- 
tion to their significance in quantum state tomography, 
SIC POMs have attracted much attention due to their 
connections with mutually unbiased bases (MUB) [THI — 
lilij , equiangular lines [2(| , Lie algebras [2l| , and founda- 
tional studies (22| . 

The trace distance is one of the most important dis- 
tance measures and distinguishability measures in quan- 
tum mechanics, and is widely used in quantum state 
tomography, quantum cryptography, and entanglement 
theory [ll. |23| - |26T | ■ as well as other contexts. In addition, 
it is closely related to other important figures of merit, 
such as the fidelity and the Shannon distinguishability 
[23|, [26[ . However, little is known about the tomographic 
resources required to achieve a given accuracy as mea- 
sured by the trace distance, since its definition involves 
taking the square root of a positive operator. Even for the 
qubit SIC POM [12, USUI]) no analytic formula is known 
concerning the mean trace distance between the estima- 
tor and the true state. One motivation of the present 
study is to solve this long-standing open problem. 

In the case of a bipartite or multipartite system, it 
is technologically much more challenging to perform full 
joint measurements such as SIC POMs on the whole sys- 
tem. Moreover, in some important realistic scenarios, 
such as tomographic quantum key distribution [29h31| , 
all parties are space separated from each other, and it 
is impractical to perform full joint measurements. Nev- 
ertheless, each party can perform a local SIC POM and 
reconstruct the global state after gathering all the data 
obtained. Such a POM will henceforth be referred to as 
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a product SIC POM; by contrast, the SIC POM of the 
whole system will be referred to as a joint SIC POM. 
The product SIC POM is particularly appealing in to- 
mographic quantum key distribution since it minimizes 
the redundant information and classical communication 
required to exchange measurement data among different 
parties (3lj . However, even less is known concerning the 
tomographic efficiency of the product SIC POM except 
for numerical studies in the two-qubit setting (2|| H2| ■ 

In this article, we aim at characterizing the tomo- 
graphic efficiency of tight IC measurements in terms 
of the mean trace distance and the mean HS dis- 
tance, with special emphasis on the minimal tight IC 
measurements — SIC POMs. We also determine the effi- 
ciency gap between the product measurements and the 
joint measurements in the bipartite and multipartite set- 
tings. 

First, we introduce random matrix theory (33| to study 
the tomographic efficiency of tight IC measurements and 
derive analytical formulae for the mean trace distance 
and the mean HS distance. We illustrate the general 
result with SIC POMs, and show the different scaling 
behaviors of the two error measures with the dimension 
of the Hilbert space. In the special case of the qubit 
SIC POM, we discuss in detail the dependence of the 
reconstruction error on the Bloch vector of the unknown 
true state and make contact with the experimental data 
given by Ling et al. [27j . As a by product, we also 
discovered a special class of tight IC measurements that 
feature exceptionally symmetric outcome statistics and 
low fluctuation over repeated experiments. 

Second, in the bipartite and multipartite settings, we 
show that the product SIC POMs are optimal among 
all product measurements in the same sense as the joint 
SIC POMs are optimal among all joint measurements. 
We further show that for bipartite systems, there is only 
a marginal efficiency advantage of the joint SIC POMs 
over the product SIC POMs. Hence, it is not worth the 
trouble to perform the joint measurements. However, 
for multipartite systems, the efficiency advantage of the 
joint SIC POMs increases exponentially with the number 
of parties. 

To provide a simple picture of the tomographic effi- 
ciency of SIC POMs and product SIC POMs, we restrict 
our attention to the scenario where the number of copies 
of the true states available is large enough to yield a rea- 
sonably good estimator, and we focus on the standard 
state reconstruction scheme, also known as linear state 
tomography [l|, |(| . The analysis of the efficiencies of other 
reconstruction schemes, such as the maximum likelihood 
method [H, [HI, [34| , is much more involved. Hopefully, our 
analysis can serve as a starting point and, in principle, 
it can be generalized to deal with those more compli- 
cated situations. Moreover, for minimal tomography on 
a large sample, the estimator given by the standard re- 
construction scheme is almost identical to that given by 
the maximum-likelihood method, except when the true 
state is very close to the boundary of the state space. 



Hence, the efficiencies of the two alternative schemes are 
quite close to each other in this scenario. 

The rest of this article is organized as follows. In 
Sec.|TTl we review the basic framework for linear state to- 
mography and recall the concepts of weighted i-designs, 
SIC POMs and tight IC measurements. In Sec. lIIH we in- 
troduce random matrix theory to study the tomographic 
efficiency of tight IC measurements, in particular the 
SIC POMs. We derive analytical formulae for the mean 
trace distance and the mean HS distance and illustrate 
the difference in their scaling behaviors with the dimen- 
sion of the Hilbert space. In Sec. lIVl we first prove the op- 
timality of the product SIC POMs among product mea- 
surements and then compare the tomographic efficiencies 
of the product SIC POMs and the joint SIC POMs. We 
conclude with a summary. 



II. SETTING THE STAGE 



A. Linear state tomography 

A generalized measurement is composed of a set of out- 
comes represented mathematically by positive operators 
Hj that sum up to the identity operator 1. Given an 
unknown true state p, the probability of obtaining the 
outcome H, is given by the Born rule: Pj — tr(ILjp). A 
measurement is IC if we can reconstruct any state ac- 
cording to the statistics of measurement results, that is 
the set of probabilities pj. If we take both the state p 
and the outcomes H, as vectors in the space of Hermi- 
tian operators, then the probability can be expressed as 
an inner product ((Hj\p)) = tr(n,p), where we have bor- 
rowed the double ket (bra) notation from Refs. [35l |36|. 
Furthermore, an out product such as ({11^ |, which is 
referred to as an superoperator henceforth, acts on this 
space just as an operator acts on the ordinary Hilbert 
space (th e arithmetics of superoperators can be found in 
Refs. [37], Hl|). With this background, one can show that 
a measurement is IC if and only if the frame superoper- 
ator 



in.MH, 
tr(n,) 



(1) 



is invertible [y, |36|, |39|, |4fj ■ The frame superoperator T 
can be written in the following form Q 



a 



where I = I !))((! I and 



inj-tr^Q/rfK^-tr^-)/^ 

tr(m) 



(2) 



(3) 



which is supported on the space of traceless Hermitian 
operators. Obviously, J- is invertible if and only if J- a is 
invertible in this space. In the rest of this article, J-q will 
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also be referred to as the frame superoperator if there is 
no confusion. 

When T is invertible, there exists a set of reconstruc- 
tion operators Qj satisfying V ■ |9j))((H,-| = I, where I 
is the identity superoperator. Given a set of reconstruc- 
tion operators, any state can be reconstructed from the 
set of probabilities pf. p — ^jPjQj- In a realistic sce- 
nario, given N copies of the unknown true state, what 
we really get in an experiment are frequencies fj rather 
than probabilities pj . The estimator based on these fre- 
quencies p — Y) j fj&j is thus different from the true 
state. Nevertheless, the deviation Ap = p — p vanishes 
in the large- N limit if the measurement is IC. In gen- 
eral, these frequencies obey a multinomial distribution 
with the MSE matrix (also called the covariance matrix) 
Sjfc = {pjSjk —PjPk) /N. The MSE matrix of the estima- 
tor p can be derived according to the principle of error 
propagation, 



c(p) = ^|e j »s ife «e fc | 

= ^(El^Mn>))(<e,Mp}}(<plV (4) 

j 

The MSE is exactly the trace of the MSE matrix, 
£u(p) = E(||Ap|| 2 HS ) = Tr(C(p)) 

^(5>tr( e i)-M? 



N 



(5) 



In this article, the symbol "Tr" is used to denote the 
trace for superoperators, and "tr" that for ordinary op- 
erators. For the convenience of later discussions, we de- 
fine N£m{p) as the scaled MSE, which is independent of 
N. The scaled mean trace distance and the scaled mean 
HS distance can be denned similarly, except that N is 
replaced by \/N. 

The set of reconstruction operators is unique for a min- 
imal IC measurement, such as a SIC POM or a product 
SIC POM, but is not unique for a generic IC measure- 
ment. Among all the candidates, the set of canonical 
reconstruction operators 



le. 



tr(n,) 



(6) 



is the best choice for linear state reconstruction in that 
it minimizes the MSE averaged over unitarily equivalent 
true states and is thus widely used in practice @. In 
the rest of this article, we will only consider canonical 
reconstruction operators. It is then straightforward to 
verify that |1)) is an eigenvector of C(p) with eigenvalue 0; 
in other words, C(p) is supported on the space of traceless 
Hermitian operators as is .7-0 . The other eigenvalues of 
C{p) determine the variances along the principle axes and 
thus the shape of the uncertainty ellipsoid. 

If TV is sufficiently large, the multinomial distribu- 
tion approximates a Gaussian distribution, which is com- 
pletely determined by the mean and the MSE matrix. 



In practice, the Gaussian approximation is already quite 
good for moderate values of N if we are mainly concerned 
with quantities like the mean HS distance and the mean 
trace distance, which are the most common figures of 
merit in quantum state tomography. We thus assume 
the validity of this approximation in the following dis- 
cussion. Under Gaussian approximation, the variance of 
the squared error |Ap|^g is given by the following simple 
formula: 



v(p) = Var(|Ap|| s ) = 2Tr(C(p) 5 



(7) 



In practice, \/v(p) quantifies the amount of fluctuation 
in the squared error |Ap|^ g over repeated experiments, 
that is the typical error in estimating £m(p) with just 
one experiment, assuming the true state is known. This 
error can be reduced by a factor of y/N e if we repeat the 
experiment iV times and take the average of ||Ap||jj S . 
In addition, once £m(p) is fixed, v(p) also quantifies the 
dispersion of the eigenvalues of C(p), that is the degree 
of anisotropy in the distribution of the estimators. 



B. Weighted i-designs and SIC POMs 



Consider a weighted set of states {\i>j),Wj} with < 
Wj < 1 and Y]j Wj — d, then the order-i frame potential 
<I> t for a positive integer t is defined as 



3,k 



S t = ^2wj{\lpj){i>j 



(8) 



Note that St is supported on the t-partite symmetric sub- 
space, whose dimension is ( 



d+t- 

t 



&t is bounded from 



and the bound is saturated if and 

l , 



below by d 2 ^ 1 - 1 

only if S t = ~ p t ym > where p t Sym is the projec- 

tor onto the t-partite symmetric subspace. The weighted 
set {\iftj},Wj} is a (complex projective) weighted t-design 
if the lower bound is saturated; it is a t-design if in addi- 
tion all the weights Wj are equal 0, Q . According to the 
definition, a weighted t-design is also a weighted i'-design 
for t' < t. 

It is known that for any positive integers d and t, there 
exists a (weighted) t-design with a finite number of ele- 
ments J4l|. However, the number is bounded from below 
by @,|4f 



d+ \t/2\ -*\(d+ [t/2\ - 1 
ft/21 ) { L*/2J 



(9) 



which is equal to d, d 2 , d 2 (d + l)/2 for t — 1, 2, 3, respec- 
tively. Any resolution of the identity consisting of pure 
states is a weighted l-design. SIC POMs @, H El and 
complete sets of MUB 0,111 El are prominent exam- 
ples of 2-designs. The complete set of MUB for d = 2 is 
also a 3-design. 
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A SIC POM is composed of d 2 subnormalized projec- 



tors onto 
fidelity 



jure states H, = \ipj){i/>j\/d with equal pairwisc 



that is, 



l^il^)| 2 = 



dS 



1 



j,k = 1,2--- ,d 2 



(10) 



It is straightforward to verify that a SIC POM is a 2- 
design from this definition. What is not so obvious is 
that a weighted 2-design consisting of d 2 elements must 
be a SIC POM j|. 

A SIC POM is group covariant if it can be generated 
from a single state — the fiducial state — under the action 
of a group consisting of unitary operators. Most known 
SIC POMs are covariant with respect to the Heisenberg- 
Weyl group (also called the generalized Pauli group) @- 
Hlj |. which is generated by the phase operator Z and the 
cyclic shift operator X defined by their actions on the 
computational basis, 



Z\e r ) = uj 
X\e r ) = 



|e r+ i) ifr = 0,l,--- ,d-2, 
|e ) ifr = d-l, 



(11) 



where ui — e 2m ^ d . A fiducial ket of the Heisenberg- 
Weyl group satisfies 



i 



Vd + T 



(12) 



for fci,fc 2 = 0,1,..., d- 1 and (ki,k 2 ) ^ (0,0). Up 
to now, analytical fiducial kets of the Heisenberg-Weyl 
CToup are known for d < 16 and d = 19, 24, 28, 35, 48 [|- 
IllL Il3| , numerical fiducial kets with high precision have 
been found up to d = 67 @, HH, 0] . 

In this article, all SIC POMs used in the numerical 
simulation are generated by the Heisenberg-Weyl group 
from the fiducial kets of Ref. |l4|. However, all theo- 
retical analysis is independent of the specific choice of 
SIC POMs. 



IC measurements when the true state is the completely 
mixed state. In that case, the MSE matrix and the MSE 
can be expressed more concisely in terms of the frame 
superoperator Tq, 



-M 



i 

N 



T 



-l 



d d 2 J 



dh^ 1 )- 



(13) 



The first equation endows the frame superoperator J-q 
with a concrete operational meaning as the inverse of the 
MSE matrix (up to a multiplicative factor) evaluated at 
the point p = 1/d. From the definitions of the frame 
superoperators F and Tq (cf. Sec. Ill A|) . we have 

Tr(Jo) = Tr(J-) - 1 < ^ tr(n,-) - 1 = d - 1, (14) 



and the inequality is saturated if and only if the measure- 
ment is rank one. Recalling that J-q is supported on the 
space of traceless Hermitian operators, whose dimension 
is d — 1, the above equation implies that 



Tr^ 1 ) > (d+l)(d 2 -l), 



(15) 



The inequalities are saturated if and only if To = Io/(d+ 
1). In other words, rank-one tight IC measurements are 
optimal in minimizing the MSE [6] . 

A rank-one tight IC measurement with outcomes H, = 
\tpj)wj(tpj\ features particularly simple canonical recon- 
struction operators 



@ j = \i> j ){d+ 



(16) 



and easy state reconstruction. According to Eq. (JSJ) , the 
MSE is also given by a simple formula [6j 



C. Tight IC measurements 

An IC measurement is tight if the frame superopera- 
tor .Fo is proportional to Io, that is J-q = alo for a > 0, 
where Io is the identity superoperator in the space of 
traceless Hermitian operators. Scott @ has shown that 
the coefficient a is upper bounded by l/(d+ 1) for any 
tight IC measurement, and the upper bound is saturated 
if and only if the tight IC measurement is rank one. He 
also showed that rank-one tight IC measurements are op- 
timal for linear state tomography in the sense that the 
MSE £m(p) averaged over unitarily equivalent density 
operators is minimized Q . Here we shall recapitulate his 
main idea in a way that suits our subsequent discussion. 

Since the average of p over unitarily equivalent states is 
the completely mixed state, according to Eqs. Q and ([5]), 
it is enough to show the optimality of the rank-one tight 



£m(p) 



N 



[d 2 +d-l- tr(p 2 )] 



(17) 



which is invariant under unitary transformations of the 
true state. In addition, the MSE matrix evaluated at 
p = 1/d is proportional to Io , which means that the un- 
certainty ellipsoid is isotropic in the space of traceless 
Hermitian operators. This feature will play a crucial role 
in our later discussions. 

There is a close relation between rank-one tight IC 
measurements and weighted 2-designs: A rank-one mea- 
surement with outcomes Uj = \^j)wj{ipj\ is tight IC if 
and only if the weighted set {\ipj}, uij} forms a weighted 
2-design [6j. For example, SIC POMs and complete sets 
of MUB are rank-one tight IC measurements according to 
this relation, which can also be verified directly. Hence, 
Eqs. (fTo]) and (TiTf are applicable to them. More exam- 
ples of tight IC measurements can be found in Ref. [7[. 
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III. APPLICATIONS OF RANDOM MATRIX 
THEORY TO QUANTUM STATE 
TOMOGRAPHY 

In this section, we apply random matrix theory to 
studying the tomographic efficiency of tight IC measure- 
ments and illustrate the general result with SIC POMs. 
In particular we derive analytical formulae for the mean 
trace distance and the mean HS distance between the es- 
timator and the true state, thus giving a simple picture 
of the resources required to achieve a given accuracy as 
quantified by either of the two distances. Our study also 
clearly shows the different scaling behaviors of the two 
error measures with the dimension of the Hilbert space. 
The idea of computing the mean trace distance using the 
random matrix theory may also be extended to derive 
other figures of merit which only depend on the devia- 
tion between the estimator and the true state. 

The rest of this section is organized as follows. In 
Sec. IIII A[ we present the simple idea of computing the 
mean trace distance and the mean HS distance with 
random matrix theory. In Sec. IIII B[ we single out 
those measurements for which the method is best jus- 
tified. In Sec. IIIICl we show that the method works very 
well for typical rank-one tight IC measurements, espe- 
cially SIC POMs. In Sec. UTTDl we focus on the qubit 
SIC POM. 



A. A simple idea 

Here is the simple idea of computing the mean trace 
distance with random matrix theory: In each experiment, 
after measurements on N copies of the unknown true 
state p, we can construct an estimator p for the true 
state according to the procedure described in Sec. Ill Al 
Once a basis is fixed, the deviation Ap = p — p can be 
represented by a d x d matrix, which varies from one 
experiment to another. After a large number of repeated 
experiments, the set of matrices Ap can be taken as an 
ensemble of random matrices obeying a multidimensional 
Gaussian distribution, which is completely determined by 
the MSE matrix C(p), 



txAp) x .■Ni,(-i((Ap|C(p)- 1 |Ap)) 



(18) 



Since C(p) is supported on the space of traceless Hermi- 
tian operators, the distribution of Ap is restricted on the 
hyperplane satisfying tr(Ap) = 0. Suppose f(x) is the 
level density function of this ensemble of matrices with 
normalization convention J dxf(x) — d. Then the mean 
trace distance between the estimator and the true state 
is proportional to the first absolute moment of f(x), 



£ tl (p) = ~E(tr|Ap|) = \ J dx\x\f(x). 



(19) 



If C(p) is (approximately) proportional to the identity 
superoperator I, then the ensemble of matrices Ap' = 



y/ d 2 /2£m(p)Ap is (approximately) a standard Gaussian 
unitary ensemble. According to random matrix theory, 
for sufficiently large d, the level density fc(x) of the 
Gaussian unitary ensemble is given by the famous Wigner 
semicircle law [331 ] : 

f (V) _ / j (2d - x 2 ) 1 / 2 if - V2d < x < V2d, (2Q) 
/gW "\ otherwise. 1 Uj 

We can derive f{x) from /g(^) by a scale transformation 
and then compute the mean trace distance between the 
estimator and the true state, with the outcome 



£tr(p) ~ —y/d£u(p)- 



3tt 



(21) 



Furthermore, one can verify that the equation is still 
quite accurate if C(p) is approximately proportional to Io 
instead of I, especially when d is large. In other words, 
the feasibility of our approach is not limited by the fact 
that C(p) is supported on the space of traceless Hermitian 
operators. 

When C(p) is proportional to Io, Ap follows a (d 2 — 1)- 
dimensional isotropic Gaussian distribution and ||Ap||^ g 
obeys \ 2 distribution with d 2 — 1 degree of freedom. The 
mean HS distance can thus be computed with the result 



£ hs (p) = E(VlAp[^ s 



£ M (p) 72r(f ) 



d? - 1 



r( 



(22) 



As a consequence of the law of large numbers, when d is 
large, £hs(p) is approximately equal to the square root 
of £m(p), and with a high probability the estimator p is 
distributed within a thin spherical shell of radius £hs(p) 
that is centered at the true state. 

In general, the accuracy of Eqs. (|2"TT) and (|2"2"j) may 
depend on the dimension of the Hilbert space and the 
degree of anisotropy of the uncertainty ellipsoid as de- 
termined by C(p). However, it turns out that the mean 
trace distance and the mean HS distance are not so sen- 
sitive to the degree of anisotropy of the uncertainty el- 
lipsoid. As we shall see shortly, the two equations are 
surprisingly accurate for a large family of measurements, 
especially tight IC measurements, even if d is very small 
(see Fig. [1]). 

Although we have started our analysis with linear state 
tomography, the idea of computing the mean trace dis- 
tance with random matrix theory has a wider applica- 
bility. We may apply the approach to study the tomo- 
graphic efficiencies of other reconstruction schemes, such 
as the maximum-likelihood method. In addition, we may 
also consider other figures of merit which only depend on 
the deviation between the estimator and the true state. 



B. Isotropic measurements 

In this section we single out those rank-one IC mea- 
surements for which the uncertainty ellipsoid is the most 
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isotropic, in which case Eqs. (j2Tj) and (1221) are best justi- 
fied. These measurements turn out to be a special class 
of tight IC measurements. In addition to minimizing the 
MSE, they also minimize the fluctuation of reconstruc- 
tion error over repeated experiments. Moreover, these 
IC measurements have the nice property that the mean 
reconstruction error is almost independent of the true 
state. 

When the true state is the completely mixed state, 
according to Sec. Ill CI C(l/d) is proportional to I if 
and only if the measurement is tight IC, and the coef- 
ficient of proportionality is minimized when the mea- 
surement is rank-one. The symmetry requirement on 
the MSE matrix is thus consistent with the efficiency 
requirement, recall that rank-one tight IC measurements 
are optimal for linear state tomography. Let us focus 
on the MSE matrix C{p) for a generic true state, assum- 
ing that we have a rank-one tight IC measurement with 
outcomes IE,- = \ipj)uij(ipj\. The degree of anisotropy 

can be quantified by Tr(C(p) 2 ) - [Tr(C(p))] 2 , where the 
over-line means taking the average over unitarily equiv- 
alent density operators. Since Tr(C(p)) is exactly the 
MSE, which is the same for all rank-one tight IC mea- 
surements according to Eq. flT7|) . it suffices to consider 
Tr(C(p) 2 ). Note that Tr (C(p) 2 ) also quantifies the fluc- 
tuation in [Ap||j S over repeated experiments according 
to Eq. 0. We find 



= d 2 



N 2 Tr(C(p) 2 

(d+l) 3 $ 3 



+ 



- 2 



2 

+ 2d - - + 
a 

2(2d 2 + 3d 



Hp 2 )] 
-i) 



2(d+l 



>d z 



2d - - 



(d-l) 
-tr(p 3 ) 4 



tr(^) 



d-l 
' d + 2 

2 r 



f 



d + 2 



(23) 



d + 2 



tr(p 2 



Hp 2 )\ 



2(d+i) 



tr(p 3 



d-l 



tr(p 2 ) 



f 



d + 2 



(24) 



d + 2 yl ' d + 2 

where $3 is the order-3 frame potential defined in Eq. ([8]). 
and we have used the inequality $3 > 6d/ (d +l)(d + 2) 
in deriving Eq. (PMl) . The lower bound is saturated if and 
only if $3 = 6d/(d+l)(d+2), that is, when the weighted 
set {\tpj),Wj} forms a weighted 3-design. 

An IC measurement derived from a weighted 3-design 
will be called an isotropic measurement for reasons that 
will become clear shortly. By virtue of the properties 
of weighted 3-designs, one can show that the MSE ma- 
trix C{p) is the same for any IC measurement derived 
from a weighted 3-design, including the covariant mea- 
surement composed of all pure states weighted by the 
Haar measure. In other words, C(p) is invariant un- 
der any unitary transformation of the measurement out- 

As an immediate consequence, 



comes: IT 



UTljW. 



the reconstruction error is the same for unitarily equiva- 
lent true states as long as the figure of merit is unitarily 
invariant, such as the mean trace distance, the mean HS 
distance, or the mean fidelity. 



Under linear state tomography, in addition to achiev- 
ing the minimal MSE £m{p), a- n isotropic measurement 
also minimizes the fluctuation of the statistical error 
over repeated experiments, or equivalently the degree of 
anisotropy in the distribution of Ap. One can show that 
with an isotropic measurement, NC(p) for a pure true 
state has only four (three if d = 2) distinct eigenvalues, 
(d+l)/(d + 2),2(d+l)/(d+2),0,2d/(d+2) with multi- 
plicities did — 2), 2(d — I), 1, 1, respectively. The degree 
of anisotropy is even lower if the true state has a lower 
purity since the leading term in the expression of C(p) 
[cf. Eq. (UJ)] is linear in p. 

In conclusion, Eqs. (f2"Tj) and (|22|) are a good approx- 
imation for computing the mean trace distance and the 
mean HS distance under isotropic measurements. After 
inserting Eq. (|17|) into the two equations, we get 



£hs(p) 



4 J d[d 2 + d - 1 - tr(p 2 ) 
3tt A 



4 d 3 / 2 



N 3ir^N 
d 2 + d - I - tr(p 2 ) V2r(f ) d 



A^(d 2 - I) 



(25) 



.(26) 



The two equations clearly show the difference in the scal- 
ing behaviors of the two error measures with the dimen- 
sion of the Hilbert space. 

An isotropic measurement is, in a sense, the most sym- 
metric measurement allowed by quantum mechanics. Re- 
markably, such a symmetric measurement can be real- 
ized with only a finite number of outcomes and its tomo- 
graphic efficiency can be characterized by simple formu- 
lae. However, since a weighted 3-design contains at least 
d 2 (d+ l)/2 elements [cf. Eq. ©], an isotropic measure- 
ment contains at least d 2 (d+ l)/2 outcomes, which are 
much more than the minimum d 2 required for an IC mea- 
surement. It is thus of more practical interests to con- 
sider generic tight IC measurements, such as SIC POMs, 
which is the focus of the next section. 



C. Tight IC POMs and SIC POMs 

In this section we consider generic rank-one tight IC 
measurements @, @] , with special emphasis on the min- 
imal tight IC measurements— SIC POMs $Ju\. When 
the weighted set Wj} forms a weighted 2-design but 

not necessarily a weighted 3-design, we can use the in- 
equality $ 3 < $ 2 = 2d/(d + I) (cf. Sec. iHBl) to derive 



an upper bound for Tv(C(p) 2 ) from Eq. (|23[) . 



N 2 Tt(C(p) 2 ) <d 2 + 2d + 2d(d + I) tr(p 2 ) - - . (27) 



In conjunction with Eqs. and (jTTJ) , this equation pro- 
vides two important pieces of information. First, the rel- 
ative deviation y/ v(p) / £m{p) is approximately inversely 
proportional to d; hence, £ns(p) is approximately equal 
to the square root of £m(p) and Eq. (|2"6")l is a good ap- 
proximation for computing the mean HS distance, espe- 
cially when d is large. Second, the degree of anisotropy 
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FIG. 1: (Color online). Theoretical and numerical simula- 
tion results on the scaled mean trace distance y/N£tr(p) and 
the scaled mean HS distance VN£hs(p) in state tomogra- 
phy with SIC POMs for d from 2 to 45. The SIC POMs 
are generated by the Heisenberg-Weyl group from the fidu- 
cial states of Ref. [l4| ]. The theoretical values are computed 
according to Eqs. (|25[) and (|26[) . respectively, with p — 1/d. 
N — 1000 + 20d 2 is chosen in the numerical simulation. The 
values for the completely mixed state is the average over 1000 
repeated experiments, and that for pure states is the aver- 
age over 1000 randomly-generated pure states each averaged 
over 100 repeated experiments. The inset shows three kinds 
of standard deviations of the scaled mean trace distances for 
the numerical simulation, (a): the standard deviation over re- 
peated experiments for the completely mixed state (the jumps 
in the curve are due to the finite number of experiments); (b): 
the average of the standard deviation over repeated experi- 
ments for each pure state; (c): the standard deviation over 
the randomly-generated pure states including a partial con- 
tribution of the fluctuation over the repeated experiments for 
each state due to the finite number of experiments. 



in the distribution of Ap cannot be too high as long as 
the measurement is rank-one tight IC. Given that the 
level density function f{x) and especially its first abso- 
lute moment are not so sensitive to slight variations in 
the degree of anisotropy, it is reasonable to expect that 
the mean trace distance can be computed approximately 
by Eq. (|25p . This expectation is supported by extensive 
numerical simulations. 

Figure Q] shows the results of theoretical calculation 
and numerical simulation on state tomography with 
SIC POMs. The mean trace distance and the mean HS 
distance from numerical simulation agree perfectly with 
the theoretical formulae in Eqs. (|25|) and (|26|) . in fact, 
much better than we have expected. The figure also 
clearly illustrates the different scaling behaviors of the 
two error measures with the dimension of the Hilbert 
space. From the inset of the figure, we see that the fluc- 
tuation in the mean trace distance over different pure 
states is much smaller than the fluctuation over repeated 
experiments on the same state. Actually, the former is so 



small that it is difficult to separate out the partial con- 
tribution of the latter with a limited number of repeated 
experiments. In other words, the reconstruction error is 
not sensitive to the identity of the true state. 

We emphasize that the results on the tomographic ef- 
ficiency of SIC POMs are representative of typical rank- 
one tight IC POMs. Since the order-3 frame potential 
$ 3 = (d 2 + 3d)/ (d + l) 2 for a SIC POM is much larger 
than the value 6d/(d + l)(d + 2) required for a 3-design, 
a SIC POM is a very poor approximation of a 3-design, 
for which Eqs. (123]) and (|26|) are best justified. Alterna- 
tively we can see this from the value of Tr(C(p) 2 ) for a 
SIC POM, which can be computed according to Eqs. ((4]), 

A 2 Tr(C(p) 2 ) = (d 2 + d + 2) [1 + tr(p 2 )] - 1 

+ [tr(p 2 )] 2 -2(d 2 + d) 2 |>3. (28) 

When d 3> 1, the term \p))((p\ in the expression of C(p) 
can be neglected and we have 



A^ 2 Tr(C(p) 2 ) « (d 2 + d) [1 + tr(p 2 )] . 



(29) 



Comparison with Eqs. (|2~4"|) and (|2~T1) shows that the value 
iV 2 Tr(C(p) 2 ) for a SIC POM is roughly in the middle 
of the lower bound and the upper bound for tight IC 
measurements. 

In the rest of this section, we briefly examine tight IC 
measurements that are not rank-one and which may arise 
in practice. In realistic experiments on quantum state 
tomography with a SIC POM, there always exists noise 
associated with detector inefficiency, dark counts, and 
other imperfections. It is important to understand how 
the noise affects tomographic efficiency. We investigate 
these effects with a simple white-noise model, in which 
the outcomes of the SIC POM are modified as follows, 



n» = 



da + d 



(30) 



where the parameter a (a > 0) characterizes the strength 
of the noise. This model is quite natural when there is 
no prior knowledge about the noise. Measurements of 
this form have also been considered in the context of 
entanglement detection with witness operators (431 ]. 

It is straightforward to verify that the measurement 
introduced above is still tight IC. The MSE can be cal- 
culated according to the procedure presented in Sec. Ill Al 
with the result 

£m( P ) = [1 + (d + l) 2 (d l)(a + l) 2 ] - tr(p 2 )}. 

(31) 

Compared with Eq. ([TT]). the MSE is roughly (a + l) 2 
times as large as in the ideal case. The mean trace dis- 
tance and the mean HS distance can still be computed 
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according to Eqs. (|2Tj) and (I221) . respectively, with the 
result 



Stv(p) 



1 



(32) 



which are roughly a + 1 times the values for the ideal 
case. Hence, due to the noise, we need roughly (a + l) 2 
times as many copies of the true states to reach the same 
accuracy as in the ideal case. A similar analysis also 
applies to tight IC measurements derived from other 2- 
designs, such as complete sets of MUB. 



D. Qubit SIC POM 



To get a concrete geometric picture, it is now better to 
work with the MSE matrix of the estimator s of the Bloch 
vector, 



C(s) 



l r 

N 



9 v 

3/ 3 ss + — {a-k ■ s) a k a k 



(34) 



where J 3 is the 3x3 identity dyadic. The mean squared 
error of the estimator s is given by 



E(|As| 2 ) = 2£ M (p) = ±(9-s 



(35) 



Suppose a 2 , a 2 , erf are the three eigenvalues of the MSE 
matrix C(s), that is the three variances along the three 
principal axes of the uncertainty ellipsoid. Then the 
mean error is determined by the following integral, 



In this section we provide further insights on the tomo- 
graphic efficiency of the qubit SIC POM by deriving an 
exact formula for the mean trace distance and discussing 
the dependence of the reconstruction error on the Bloch 
vector of the true state (see Refs. [HI, H3, [HI for earlier 
accounts) . We also confirm that the result based on ran- 
dom matrix theory is already quite accurate for d = 2, 
although it is best justified when d is large. As a sim- 
ple application, we make contact with the experimental 
result given by Ling et al. [27j |. 

For the qubit SIC POM, the four outcomes U k for k = 
1,2,3,4 are in one-to-one correspondence with the four 
unit vectors a k pointing to the four vertices of a regular 
tetrahedron inscribed on the Bloch sphere; that is, = 
j(l + a k -t), where t — (r 1; r 2 , r 3 ) are the Pauli matrices 
(<7j are reserved to denote the standard deviations in this 
article). The reconstruction operators are given by Q k = 
5(1 + 3<2fc • t) according to Eq. ([16]). Let s denote the 
Bloch vector of the true state p. To reconstruct the true 
state p is equivalent to reconstruct the Bloch vector s 

a 

4 1/ 4 \ 



k=l 



k=l 



3}p k a k . 

fe=i 



(33) 



Meanwhile, both the HS norm |A/j[hs and the trace 
norm |Ap|t r are proportional to the Euclidean length 
of As = s — s, where s is an estimator of s; that is, 
lApllg = (As) 2 /2, ||Ap|| tr = \As\/2. 

The MSE matrix of the estimator p can be calculated 
according to Eq.(jU), with the result 



C{ P ) 



3=1 



S ■ T S • T 



16 



53 \ a k ■ T )) a k ■ s{{a k ■ t\ 



EflAsI 



dxdydz 



\J x 1 + y 2 + z 2 
(27r) 3 / 2 (Ticr2cr3 



x exp 



2o\ 



2a 2 



2a 2 



3 3/2 



where 

9 = ^-t 2 ) 2 



'l u 2 .4 



< 4 + K+a 2 )t 2 (l-t 2 ). 



If at least two of the variances are equal, say a<i = 01, 
then the integral can be evaluated explicitly, 



E(|A fl |) 



-cr 3 

7T 



2 ( a \ arctan y 
2" /tr 2 arctanh 



if cti = 0, 
ifa 3 =0, 
if CT 3 = (Ti, 

0-3) ifci > CT 3 , 
cr 3 I if (Ti < cr 3 . 



(36) 



If the uncertainty ellipsoid is isotropic, that is, a± 
&2 = c 3 , then we have 



Strip) 



(37) 



fe=i 



For the completely mixed state, this equation is exact, 
by contrast, the expression £tr(p) ^ 4-\/9 — s 2 / (Bny/N) 
based on random matrix theory [see Eq. pip ] is about 8% 
smaller. The disparity is much smaller than the relative 
deviation of |Ap|t r over repeated experiments, which is 
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FIG. 2: The scaled mean trace distances for states with Bloch 
vectors that are either parallel or antiparallel to the legs of the 
qubit SIC POM (theory), and the scaled mean trace distance 
averaged over randomly-generated states with a fixed purity 
(numerical simulation). N = 1000 is chosen in the numeri- 
cal simulation, and the scaled trace distance for each given 
purity is averaged over 1000 randomly-generated states each 
averaged over 400 repeated experiments. The dotted curve 
shows the scaled mean trace distance determined by Eq. (|37[1 , 
which assumes that the uncertainty ellipsoid is isotropic. 



states sits roughly in the middle of the two extreme cases. 
In all three cases, there is a slight decrease in the mean 
trace distance as the purity of the true state increases, 
which can roughly be attributed to two reasons: the de- 
crease in the MSEs [cf. Eq. (f35j) ] and the increase in the 
degrees of anisotropy of the uncertainty ellipsoids. 

Ling et al. 27] have studied the tomographic effi- 
ciency of the qubit SIC POM experimentally and de- 
termined the scaled mean trace distances for the three 
states with z = 0,-1,1, respectively, with the result 
1.417,1.288,1.323. By contrast, our theoretical calcu- 
lation has yielded the result 1.382, 1.259, 1.295. The ex- 
perimental and the theoretical values reflect the same de- 
pendence of the reconstruction error on the Bloch vector 
of the true state. The former are slightly larger than the 
latter, but the difference is very small, in other words, 
the agreement between experimental data and theoret- 
ical calculation is pretty good. Note that the relative 
fluctuation of the reconstruction error over repeated ex- 
periments is larger than 40% and the experimental values 
are the average of only 40 runs. In addition, any imper- 
fection inevitable in real experiments may also affect the 
accuracy of the estimator. 



about 42%. For other true states, the disparity is even 
smaller. Hence, the result based on random matrix the- 
ory is already quite accurate even for d = 2. 

Those states whose Bloch vectors are either parallel or 
antiparallel to the legs of the SIC POM have attracted 
more attention both theoretically [l2j and experimentally 
[27) . since they represent two extreme cases. We shall 
compute the mean trace distances for those states and 
discuss the dependence of the reconstruction error on the 
Bloch vector of the true state. 

If ai is chosen as the z axis, without loss of generality, 
the Bloch vectors of those extreme states can be param- 
eterized as s = za\ with — 1 < z < 1. According to 
Eq. (|3"4")l . the MSE matrix of s now reads 



C(s) 



(3 — z)I 3 + (3z — z 2 )a,ia,i 



N 



whose eigenvalues are given by 



N 



(3 -*)(! + *) 
N 



(38) 



(39) 



Note that the uncertainty ellipsoid is rotationally sym- 
metric. As z decreases from 1 to —1, the uncertainty 
ellipsoid evolves from prolate to oblate and finally to a 
singular ellipsoid when z = — 1. The mean trace distance 
of those extreme states can be calculated according to 
Eq. (|36[), Figure [5] shows the scaled mean trace distances 
for those states together with the numerical average over 
randomly-generated states. The mean trace distance is 
slightly smaller for states with Bloch vectors that are an- 
tiparallel to the legs of the SIC POM than those that are 
parallel. For a fixed purity of the true states, the aver- 
age of the mean trace distance over randomly-generated 



IV. JOINT SIC POMS AND PRODUCT 
SIC POMS 

In the bipartite or multipartite settings, it is tech- 
nologically much more challenging and sometimes even 
impossible to perform full joint measurements such as 
SIC POMs on the whole system. It is thus of paramount 
practical interests to determine the optimal product mea- 
surements and the efficiency gap between the product 
measurements and the joint measurements. We show 
that under linear state tomography, product SIC POMs 
are optimal among all product measurements in the same 
sense as joint SIC POMs are optimal among all joint 
measurements. Furthermore, in the bipartite setting, 
there is only a marginal efficiency advantage of the joint 
SIC POMs over the product SIC POMs and it is thus 
not worth the trouble to perform the joint measure- 
ments. However, in the multipartite settings, the effi- 
ciency advantage of the joint SIC POMs over the prod- 
uct SIC POMs increases exponentially with the number 
of parties. 



A. Bipartite SIC POMs and product SIC POMs 

Consider a product measurement on a bipartite sys- 
tem whose parts A, B have subsystem dimensions d\ , c?2 
respectively, and the total dimension is d = d\di- To 
show the optimality of the product SIC POM, we shall 
use the same strategy described in Sec. Ill Cl More gener- 
ally, we show that if the product measurement minimizes 
the MSE averaged over unitarily equivalent states, then 
the measurement on each subsystem is rank-one tight 
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IC, and vice versa. As an immediate consequence, the 
product SIC POM is optimal and furthermore any min- 
imal optimal product measurement must be a product 
SIC POM, recall that SIC POMs are the only minimal 
tight IC measurements [6(. 

Since the average of p is the completely mixed state, it 
suffices to demonstrate our claim when p = 1/d according 
to Eq. ((4]) . Suppose IL^ are the outcomes of the measure- 
ment on the first subsystem and IL- 2 are those for the sec- 
ond subsystem, then each outcome in the product mea- 
surement has a tensor product form Tlj 1 j 2 = ® n j2 . 
The same is true for the frame superoperator J- = J-\®!F2 
and the reconstruction operators Qj 1 j 2 = Qj 1 <8> Qj 2 - Ac- 
cording to Eq. (fl"3"|). we have 

/In _ Jx 

£m Q = ^[ Tt (- F r 1 ) Tr (^" 1 )-i]- (40) 

The MSE is minimized if and only if both T^J 7 ^ 1 ) and 

Tr^J 7 ^ 1 ) are minimized, that is, when the measurement 
on each subsystem is rank-one tight IC (cf. Sec. Ill C|) . 

Next, let us focus on the tomographic efficiency of 
the optimal product measurements. If the product mea- 
surement is composed of two rank-one tight IC measure- 
ments, as in the case of the product SIC POM, then each 
factor in the reconstruction operator JU2 = O^ <B>0j 2 is 
given by Eq. (Q~6|) . The MSE can be computed according 
to Eq. ©, 

^° d (p) = ^ M + di - l)(dl + da - 1) - t V)] . (41) 

Surprisingly, the MSE is almost independent of the true 
state, as in the case of SIC POMs. In addition, it is ap- 
proximately equal to the product of the MSEs for two 
subsystems, respectively. The variance v prod (p) of the 
squared error may depend on the specific choice of prod- 
uct measurements according to Eq. ([7]). For the product 
SIC POM, it is approximately given by 

vP1 '° d (p) * ^2 {( d i + d i - 2 )( d 2 + d 2 - 2) 

x [1 + tr(p 2 ) + tr(pj) + tr(pl)] + (dj + d 1 - 2) 
x [1 + tr(p 2 )] + (dj + d 2 - 2) [1 + tr(p 2 )] }. (42) 

Note that the variance does not only depend on the pu- 
rity of the global state, but also on the purities of the 
reduced states, which means that it generally depends 
on the entanglement of the global state. For example, if 
the true state is pure, the variance is approximately max- 
imized for product states and minimized for maximally 
entangled states. 

Compared with the result of the joint SIC POM given 
in Eq. JTZJI, the MSE achieved by the product SIC POM 
is slightly larger, but the difference is generally very 
small, especially when both d\ and di are large. By 
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FIG. 3: The ratio of the MSEs in state tomography with the 
product SIC POM and the joint SIC POM when the true 
state is the completely mixed state, note that the ratios for 
other true states are almost the same. The ratio is maximized 
when di = d<2 = 3. The ratio of the mean trace distances is 
approximately equal to the square root of the ratio of the 
MSEs. 

contrast, the fluctuation over repeated experiments is 
stronger by a bigger margin in the product SIC POM. 
Figure [3] shows the ratio of the MSEs when the true state 
is the completely mixed state, note that the ratio for 
other true states is almost the same. The maximal ratio 
1.36 is obtained when d± = di = 3. If d\,di > 3, the ra- 
tio decreases monotonically with d\ and dn ; if d^ = 2 and 
d\ > 3, the ratio decreases monotonically with d\. For 
sufficiently large di,d%, the ratio is about l + l/di + l/ cfo. 
In conclusion, there is only a marginal efficiency advan- 
tage of the joint SIC POM over the product SIC POM. 
The product SIC POM is thus more appealing for practi- 
cal applications since it is much easier to be implemented. 

Although the product SIC POM is not even a tight IC 
measurement, comparison of Eqs. (|4Tj) and (|42|) shows 
that the relative deviation of the squared error is quite 
small, especially when di,c?2 are large. Hence, Eq. (|2"2")l 
is still a good approximation for computing the mean 
HS distance. The mean trace distance can be calculated 
approximately according to Eq. (|21[) . with the result 

Strip) * \A d ? + dl - ^ + d 2 - 1) - tr(p 2 ). 

(43) 

Generally speaking, the larger the values of d\ and di 
are, the more accurate is this formula. The ratio of the 
mean trace distance with the product SIC POM to that 
with the joint SIC POM is thus approximately equal to 
the square root of the ratio of the MSEs. 

Table |l] shows the theoretical and numerical simula- 
tion results on the scaled mean trace distances for the 
two-qubit product SIC POM and joint SIC POM. There 
is quite a good agreement between theoretical calcula- 
tion and numerical simulation although eq and di are so 
small. The mean trace distances achieved by the product 
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TABLE I: Theoretical and numerical simulation results of the 
scaled mean trace distances for the two-qubit joint SIC POM 
(Joint) and product SIC POM (Prod). The theoretical values 
are computed according to Eqs. (|25[1 and (|43[) . respectively. 
N = 1000 is chosen for the numerical simulation. For the 
completely mixed state, the scaled trace distance is averaged 
over 1000 repeated experiments. For pure states, it is av- 
eraged over 1000 randomly-generated states, each averaged 
over 1000 repeated experiments. The standard deviations of 
the scaled trace distances over the 1000 randomly-generated 
pure states (including a partial contribution of the fluctuation 
over the repeated experiments for each state due to the finite 
number of experiments) are 0.033 and 0.027 for the product 
SIC POM and the joint SIC POM, respectively, both of which 
are very small. 



POM 


Completely mixed state 
Theory Numerical Error% 


Average over pure states 
Theory Numerical Error% 


Prod 
Joint 
Ratio 


4.223 4.255 -0.8 
3.676 3.716 -1.1 
1.149 1.145 


4.158 4.162 -0.1 
3.601 3.575 +0.7 
1.155 1.164 



SIC POM are roughly 15% larger than that achieved by 
the joint SIC POM. As a consequence, with the product 
SIC POM, we need about 32% more copies of the true 
states to reach the same accuracy achieved by the joint 
SIC POM. Despite its slightly lower efficiency, the prod- 
uct SIC POM is more appealing than the joint SIC POM 
due to the relative ease in its implementation in real ex- 
periments. The same conclusion has also been reached 
in Ref. (32[, where the maximum-likelihood method was 
adopted for state reconstruction. 



B. Multipartite SIC POMs and product SIC POMs 

Suppose k parties want to reconstruct a quantum state 
shared among them with a product measurement, and dj 
for j = 1, 2, • • • , k is the dimension of the Hilbert space 
of the jth party. According to the same analysis as in 
the bipartite setting, under linear state tomography, the 
product SIC POM is optimal among all product measure- 
ments. The MSE achieved by the product SIC POM can 
also be calculated in the same manner, with the result 



~prod 

'm 



(p) 



1 

N 



■ k 

JJ^ + df-lJ-tr^) 

■3=1 



(44) 



which is almost independent of the true state. When the 
true state is the completely mixed state, the variance of 
the squared error is given by 



k 

n 

3=1 



(d? + 2d] - 2) 



k 

3 = 1 3 ■ 



(45) 



For a generic true state, the variance depends on the 
purity of the global state as well as the purities of various 
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FIG. 4: (Color online). Theoretical and numerical simula- 
tion results of the scaled mean trace distances for the joint 
SIC POMs and the product SIC POMs on multiqubit sys- 
tems, where k is the number of qubits. The theoretical values 
are computed according to Eqs. (|25|) and (|47[) , respectively, 
with p = 1/d. N — 1000 + 20d 2 is chosen in the numeri- 
cal simulation. For the completely mixed state, the scaled 
trace distance is averaged over 1000 repeated experiments. 
For pure states, it is averaged over 1000 randomly-generated 
states, each averaged over 100 repeated experiments. 



reduced states, and can be much larger than the value 
given above. 

If the dimension of the Hilbert space of each party is 
equal to g?i , then the ratio of the MSE for the product 
SIC POM to that for the joint SIC POM grows exponen- 
tially with the number of parties k, 



~prod 
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(p) 



p joint 



(P) 
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1 

d~i 



1 

dj 



(46) 



The ratio of the variances grows with an even higher 
rate, whose specific value may heavily depend on the 
true state. In other words, the efficiency advantage of 
the joint SIC POM over the product SIC POM grows 
exponentially with the number of parties. 

Although the fluctuation in the reconstruction error 
over repeated experiments is stronger in the product 
SIC POMs as compared with the joint SIC POMs, the 
relative fluctuation is still small. Hence, Eq. (|22t is still 
a good approximation for computing the mean HS dis- 
tance. When k is not too large, the mean trace distance 
can be calculated approximately according to Eq. (|21l) . 
with the result 



£trO) 



AVd 



3tt 



m 

3 = 1 



d,-l)-tr(p 2 ). (47) 



Since the ratio of the mean trace distance achieved by the 
product SIC POM to that achieved by the joint SIC POM 
is approximately equal to the square root of the ratio of 
the MSEs, it also increases exponentially with the num- 
ber of parties; the same is true for the mean HS distance. 
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Figure [4] shows theoretical and numerical simulation re- 
sults of the scaled mean trace distances for the product 
SIC POMs and the joint SIC POMs on multiqubit sys- 
tems (numerical fiducial kets of SIC POMs are available 
at Ref. [14J). There is a pretty good agreement between 
theoretical prediction and numerical simulation for k up 
to 5. This plot further confirms that the efficiency advan- 
tage of the joint SIC POM over the product SIC POM 
increases exponentially with the number of parties. 

V. SUMMARY 

We have introduced random matrix theory [331 ] for 
studying the tomographic efficiency of tight IC measure- 
ments, which include SIC POMs as a special example. In 
particular, we derived analytical formulae for the mean 
trace distance and the mean HS distance between the 
estimator and the true state and showed the different 
scaling behaviors of the two error measures with the di- 
mension of the Hilbert space. The accuracy of these 
formulae were confirmed by extensive numerical simu- 
lations on state tomography with SIC POMs. In the 
special case of the qubit SIC POM, we derived an ex- 
act formula for the mean trace distance and discussed in 
detail the dependence of the reconstruction error on the 
Bloch vector of the unknown true state. As a byprod- 
uct, we also discovered a special class of tight IC mea- 
surements called isotropic measurements, which feature 
exceptionally symmetric outcome statistics and low fluc- 
tuation over repeated experiments. 

In the bipartite and multipartite settings, we showed 



that the product SIC POMs are optimal among all 
product measurements in the same sense as the joint 
SIC POMs are optimal among all joint measurements. 
We further showed that for bipartite systems, there 
is only a marginal efficiency advantage of the joint 
SIC POMs over the product SIC POMs, which disap- 
pears in the large-dimension limit. Hence, it is not worth 
the trouble to perform the joint measurements at the 
current stage. However, for multipartite systems, the ef- 
ficiency advantage of the joint SIC POMs over the prod- 
uct SIC POMs increases exponentially with the number 
of parties. 

Our study provided a simple picture of the scaling 
behavior of the resource requirement in state tomogra- 
phy with the dimension of the Hilbert space, and of the 
efficiency gap between product measurements and joint 
measurements. The idea of applying random matrix the- 
ory to studying tomographic efficiencies may also find 
wider applications in other state estimation problems. 
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